Vortices in the presence of a nonmagnetic atom impurity in 2D XY ferromagnets 
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Using a model of nonmagnetic impurity potential, we have examined the behavior of planar 
vortex solutions in the classical two-dimensional XY ferromagnets in the presence of a spin vacancy 
localized out of the vortex core. Our results show that a spinless atom impurity gives rise to an 
effective potential that repels the vortex structure. 
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The easy-plane Heisenberg ferromagnet in two dimen- 
sions and continuum limit supports non-linear pseudo- 
particles with a vortex structure. These excitations are 
of paramount importance in the understanding of static 
and dynamical properties of magnetism. For example, 
the vortex unbinding is responsible for a phase tran- 
sition known as Kosterlitz-Thouless transition Be- 
sides, these may be responsible for a central peak in the 
dynamical correlation function ^ |^, ||, || observed in 
Monte-Carlo simulations ^, |?| and experiments ||. The 
simplest vortex configuration, referred as the planar vor- 
tex, occurs when the anisotropy is large, resulting in spin 
confinement to the lattice plane || |l0) . 

The interaction of vortices with spatial inhomo- 
geneities is of considerable importance from both the 
purely theoretical and applied point of view. Impurities 
and/or defects are present even in the purest of material 
samples and their effect on the motion or structure of 
nonlinear excitations must be considered when the dy- 
namics or configurations of such solutions are important 
in the problem at hand. Recently, Zaspel, McKennan 
and Snaric investigated, using the discrete lattice, 
the instability of planar vortices and concluded that they 
will be stable at a large range of anisotropics if there is 
a nonmagnetic impurity such as Cd or Zn at the center 
of the vortex. In this paper we study, using the con- 
tinuum approximation, the interaction between a planar 
vortex and a nonmagnetic impurity localized out of the 
vortex center. To this end, we start defining the classical 
XY ferromagnet model, which is given by the following 
Hamiltonian: 



H 



(1) 



Where J is a coupling constant, the classical spin vec- 
tor has three components S = (S x , S y , S z ) and the sum- 



mation is taking over the nearest-neighbor square lattice 
sites. This model is one of the most studied in statis- 
tical physics and has been found to describe a wide va- 
riety of systems with complex scalar order parameters, 
including superconductor films, Josephson junction ar- 
rays, and superfluid He 4 films. The choice of the XY 
model is arbitrary for our purpose, since the results 
can be used, without modifications, to any other model 
with XY symmetry such as the classical easy-plane fer- 
romagnetic model. It is convenient to parametrize the 
spin field in terms of spherical coordinates as follows: 
S = S(cos 8 cos <fi, cos 6* sin </>, sin 9). By a straightfor- 
ward generalization of arguments used to obtain the con- 
tinuum limit of Heisenberg Hamiltonians in the case of 
one dimension we can write the continuum version 
of the Hamiltonian ([!]) as 



d 2 r 



m 2 (Vm) 2 



(1 - m 2 )(V$) 2 + -^m 2 



(2) 
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where m — sin 6, ao is the lattice constant and we 
have taken S 2 = 1. One can obtain the motion equa- 
tions, dm/dt = SH/5<f>, d<S>/dt = -SH/5m, for this 
field theory in the usual way using the pair of canoni- 
cally conjugated variables m and 

Now, consider that the system contains a nonmagnetic 
impurity concentration in the plane. For simplicity, we 
consider only one nonmagnetic atom present. If we re- 
move a spin from the lattice, the nearest neighbors of that 
spin will have a coordination number of three, instead of 
bulk-spin coordination number of four. Therefore, such 
boundary spins would have larger fluctuations than the 
bulk spins and it is conceivable that the non-linear config- 
urations such as a vortex, would preferentially nucleate 
around this vacancy. In fact, in this circumstance, the 
vortex energy is lowered, since the nonmagnetic impu- 
rity at the vortex center will remove the nearest-neighbor 
exchange bonds at the impurity in a radially symmetric 
way without modifying the symmetric configuration of 
the vortex, while its energy in the region without the 
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impurity remains the same. However, the vortex energy 
increases logarithmically with the system size L, and in 
an infinitely extended system, this energy would diverges 
as L — > oo so that we should not expect that a single 
vortex could nucleate around the spin vacancy. In fact, 
vortices are always created in pairs of vortex-antivortex 
having finite energy and separated by a few lattice con- 
stants. Then, if one member of the pair nucleates around 
the impurity, the other member will be near it. Since a 
vortex pair does not have a cylindrical symmetry, the 
energy of this system does not necessarily decrease, al- 
though there are less nearest-neighbor exchange bounds, 
because the spin vacancy may deform the pair configura- 
tion, increasing its energy. Nevertheless, the interest of 
this paper is to study the behavior of a vortex in the sys- 
tem, that is not nucleated at the position of the vacancy. 
In this case things must change considerably because the 
spin vacancy may also deform the single vortex config- 
uration. To take into account the nonmagnetic impuri- 
ties we consider the following modified XY ferromagnetic 
Hamiltonian in the continuum limit 



/// L 2 l 



m 2 (Vm) 2 

1 — TO 2 



(1 -m 2 )(V$) 2 + -^m 2 



by 



xV(f),(3) 

where V(f) is a nonmagnetic impurity potential given 



One point to note in this equation is its dependence 
on the spin field around the position of the impurity. 
If the vacancy is localized in a region where the spin 
configuration consists of aligned spins like a domain with 
all spins aligned along some direction (|V$| « 0), the 
spin field practically does not feel the presence of the 
impurity. However, an impurity placed in a region where 
the spin directions vary considerably (|V$| >> 1/ao) 
may have a strong coupling with the spin field and may 
modify the initial spin configuration for large distances. 

In polar coordinates the vectors f and fo are written 
as (r, 0), and {r ,4>o) respectively. To solve Eq.(7), we 
note first that the gradient of the impurity potential can 
be expressed as 



W(f) = a [rcos(a - \<p - <j) Q \) + 0sin(a -\<j>- <j> \)} 

xS(f-f -b),(8) 

where d is the Dirac delta function and a is the angle 
that the vector b, with origin at the point fo and end at a 
point on the circumference of the potential, makes with 
the vector fo. As we are interested in a local impurity 
with atomic dimensions, we make b — > in the contin- 
uum limit (to be more precise, we should make b — > do) 
indicating that the impurity is an atom (such as Zn, Mg 
or Cd for example). In this case, f — > fo, 4> — » 4>o, and 
we rewrite Eq. (0) as 



V(r) 



1 if \f — fo\ >b, 
if f - f < b. 



(4) 



Here the impurity is centralized at the point fo and has 
the form of a circle with diameter equal 26. There is a 
circular region in the plane, around the point fo, without 
any magnetic interaction. 

Substituting Eq. (|]) into the equations of motion we 
get 



1 d6 

-— = cos 8V(f)V 2 <S> - 2 sin 8V(f)V9 • V$ 



-cos 0W(f) ■ V$, (5) 



W(f) ~ ao[rcos(a) + 4>sm(a)]5(f — fo), 



(9) 



where we can interpret cos(a) and sin(a) as anisotropic 
coupling constants. This coupling depends on the direc- 
tion one looks, if the observer center is placed on the 
impurity position. 

Considering Eq. ([?]) with V(f) = 1 at the left side 
(this fails only at the point fo, since the impurity is local) 
and supposing that the vortex structure is modified by 
the presence of the nonmagnetic impurity, we write $ = 
$0 + ^1, where $0 = arctan(y/x) is the traditional single 
vortex solution for a vortex with center localized at the 
origin and $1 is the deformation caused by the spinlcss 
impurity localized at fo. Thus, Eq. ^ with the above 
considerations can be written as 



1^ = - tan 6 sin 6V{f)V 2 6 - sin eV{f){V6) 2 

+ sin(9V7f)[4/a 2 - (V$) 2 ] - tan 2 9 sin 6>W(f) • V6. (6) 

Our interest is in planar and static solutions in the 
presence of this nonmagnetic impurity potential. Hence, 
we take dO/dt = d&/dt = and m = sin 8 = in Eqs. 
(H) and (H), obtaining only one and simpler equation to 
be solved 

T/(f)V 2 $ = -W(f) • V$. (7) 



V 2 ($ + $1) = -a V($ + $1) 

■[f cos(a) + 4>sin(a)]5(f — fo) (10) 

Using the fact that V 2 $ = and taking V^o + '&i) = 
V'I'o = f$ near the point fo, Eq. ( |l0| ) can then be ap- 
proximated by 

V 2 $! = -^sinfaWf -f ), (11) 

or 
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-27rro<i>i 



clq sin(a) 



2ir5(r — r*o). 



(12) 



This is easily solved using the fact that in two dimen- 
sions, V 2 ln(r) = 2ir5(r). We get 



ao 



(13) 



Writing the anisotropic coupling constant along the a- 
direction in terms of r, <p, the vortex structure with its 
center at the origin in the presence of a nonmagnetic 
impurity localized at r*o is given by 



<& = arctan(?//a;) 



ao r sm( 



») 



2nr 



In 



r 



\r - r \ 
a 



(14) 



The configuration of this deformed vortex is shown in 
figures 1 and 2. Although the continuum theory can- 
not be applied near the vortex core, in Fig. 1 we have 
considered the impurity one lattice sapcing from the vor- 
tex center just to emphasize the vortex deformation as 
the vortex core approaches the impurityWe notice that 
if ro is large (the vortex center is far away from the spin 
vacancy) the vortex practically keeps the same original 
form, but for small ro the vortex configuration suffers a 
severe modification, mainly in the region in which the im- 
purity is located. This is due to the fact that the gradient 
of the spin field is small in the region of the impurity if it 
is far way from the vortex center and large in the region 
of the impurity if it is near the vortex center. When an- 
impurity is near the vortex core, Eq. ( pT| ) implies (as it 
can also be seen partially in figure 1) that a large domain 
with all spins aligned along the direction perpendicular 
to fo will be formed in a region located after the impurity 
position ro- 
In order to calculate the energy of this planar solution 
we consider the Hamiltonian (|^) with m — 0, obtaining 
Ei = f(VQ) 2 V(r)d 2 r. As we have seen, the field $ 
describes a single vortex at the origin in the presence of 
an impurity at distance ro away. The effective potential 
experienced between the two defects (one defect in the 
spin field and the other in the lattice structure) is defined 



U eff (r )=E I -E 1/ , 



(15) 



where E v = 7rJln(L/ao) is the energy of a single vor- 
tex in the absence of impurities. Making suitable approx- 
imations, we find that such effective potential results in a 
repulsive central interaction with a dominant term given 
by 



a%El 



1 



24tt 4 J 2 r 2 



(16) 



We see, therefore, that the presence of a nonmag- 
netic impurity increases the vortex energy as the dis- 
tance between the impurity and the vortex decreases. 
In a ferromagnet with a size of the order of L sw 10 8 ao 
( a few centimeters), a spinless atom impurity situated 
about 2ao from the vortex core would increase the vor- 
tex energy about 36%. Note that the effective potential 
barrier becomes infinity as ro — > and it is energeti- 
cally favorable that vortices and impurities become far 
apart. But, if the calculations were taken considering 
that the spin vacancy is localized at the vortex center, 
we would have W(r) = ao5(f)r and near the vortex 
core V$ « (l/a )(f>, leading to y(r)V 2 $ = 0. As a 
consequence, in the region without the spinless impu- 
rity, where V 2 <& = 0, one gets the same typical solutions 
and the vortex structure does not suffer any alteration. 
Hence, the only one effect of a central nonmagnetic im- 
purity, is to make the vortex energy to decrease, because 
of the nonexitence of nearest-neighbor exchange bonds 
at the impurity. Nevertheless, as we suggested earlier, a 
single vortex with infinite energy may not nucleate by it- 
self around the impurity, since these are created in pairs. 
Besides, Eq. ( |l6| ) shows that an infinite potential bar- 
rier has to be exceeded by the vortex core in order that 
it might reach the nonmagnetic impurity and the mini- 
mum of energy. Then, one should not expect to find a 
single vortex with a spin vacancy localized in its center. 

In Summary, vortices prefer to stay far way from non- 
magnetic impurities and hence, the spin dynamics must 
be affected by these lattice defects. Our calculations 
could also be taken for two-dimensional easy-plane an- 
tiferromagnets. It would be carried out in essentially 
the same way, leading to similar results. The structure 
and motion of vortices in two-dimensional magnets may 
be driven by the presence of spinless impurities due to 
the repulsive effective potential. Since the dynamical 
structure factor is the Fourier transform of the vortex 
spatial and temporal configuration, we expect that non- 
magnetic impurities may cause changes in the central 
peak j| H |, fl and also in the electron paramagnetic 
resonance linewidth ]l3[ ], which must be seen in neutron 
scattering and resonance experiments. However, much 
work has to be done in order to see these effects. Be- 
sides, since the vortex energy is modified, the Kosterlitz- 
Thoulcss temperature may also be affected by the pres- 
ence of impurities. In fact this theory holds also for 
temperatures below the Kosterlitz-Thouless temperature 
Tkt, where vortices are bound in pairs. However, the 
problem of vortex pairs interacting with nonmagnetic 
impurities and its influences on Tkt will be treated in 
a future paper.We also suggest that the above calcula- 
tions may have some relevance to high-T c superconduc- 
tors, because a common feature of all high-transition- 
temperature cuprates is the proximity between antiferro- 
magnetic and d-wave superconducting phases controlled 
by the doping. The effect of impurities on supercon- 
ductors has been of theoretical and experimental interest 
even in its own right for a long time. Recent nuclear mag- 
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netic resonance measurements have shown that when a Acknowledgments 
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FIG. 1: Structure of a single vortex with center at (0,0) in the presence of a nonmagnetic impurity located at the site (1,0). 
Since the nonmagnetic impurity is near the vortex center, the vortex experiences a strong effect of the spinless atom impurity. 
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FIG. 2: Here, the impurity is located at (5,0). Note that the vortex structure is almost perfect, indicating that a vacancy put 
away from the vortex core has small influence on the vortex structure. 



